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^> ' I briefly review the results of recent searches for flavor-changing neutral current and 

QQ ^ lepton-flavor and lepton-number violating decays of D"*", -D3 , and D" mesons (and their 

antiparticles) into modes containing muons and electrons. The primary focus is the 

results from Fermilab charm hadroproduction experiment E791. E791 examined 24 7t££ 

^ ' and KU decay modes of D"*" and D^ and £'^£~ decay modes of D^. Limits presented by 

r-^ ^ E791 for 22 rare and forbidden dilepttm decays of D mesons were more stringent than 

^vj ■ those obtained from previous searches,E| or else were the first reported. 
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^P ' 1. Introduction 

O , 

^5 , Searches for rare and forbidden decays allow investigation of physics beyond the 

QJ ' Standard Model and of phenomena in mass ranges beyond those available to cur- 



rent accelerators. One way to search for physics beyond the Standard Model is to 
(D I examine decay modes that are forbidden or else are predicted to occur at a negligi- 

i-^ ' ble level (rare). Observing such decays would constitute evidence for new physics, 

and measuring their branching fractions would provide insight into how to modify 
our theoretical understanding, e.g., by introducing new particles or new gauge cou- 
}J] ' plings. These new mediators could include: leptoquarks, horizontal gauge bosons, 

5^ , or something even more esoteric. 

This review focuses on recent results from the Fermilab experiment E791a search 
for 24 decay modes of the neutral and charged D mesons (which contain the heavy 
charm quark). These decay modes* fall into three categories: 

(i) FCNC - flavor-changing neutral current decays {D^ -^ £+£^ and Dt^ s -^ 

h'^l^l^, in which h is either a tt or K); 
(ii) LFV - lepton-flavor violating decays (Z?'' -^ /-t^e^, Df^ ^. -^ /i^/i^e^, and 
DJ,^ ^^ -^ /i^/i+e+, in which the leptons belong to different generations); 

* Charge-conjugate modes are included implicitly throughout this paper. 



(iii) LNV - Icpton- number violating decays {D'\ ^^ -^ h £+^+, in which the lep- 
tons belong to the same generation but have the same sign charge). 

Decay modes belonging to (i) occur within the Standard Model via higher-order 
diagrams, but estimated branching fractionsd are 10~^ to 10^^. Such small rates 
are below the sensitivity of current experiments. However, if additional particles 
such as supersymmetric squarks or charginos exist, they could contribute additional 
amplitudes that would make these modes observable. Decay modes belonging to 
(ii) and (iii) do not conserve lepton number and thus are forbidden within the 
Standard Model. However, lepton number conservation is not required by Lorentz 
invariance or gauge invariance, and a numbeiiof theoretical extensions to the Stan- 
dard Model predict lepton-number viQlation.a Many experiments have searched for 
lepton-number violation in K dccaviLLl and fonJepton-number violation and flavor- 
changing neutral currents in D at3 and B t3 decays. While many experiments 
have examined charge 1/3 strange and beauty quarks, investigations of D mesons 
look for rare and forbidden decays involving charge 2/3 charm quarks. Charge 2/3 
quarks may couple differently than charge 1/3 quarksD 

2. E791 Analysis 

The data were gathered by the the E791 spectrometer Jl3 which recorded 2 x 10^° 
events with a loose transverse energy trigger. These events were produced by a 500 
GeV/c TT^ beam interacting in a fixed target consisting of five thin, well separated 
foils. Track and vertex reconstruction was provided by 23 silicon microstrip planeaiil 
and 45 wire chamber planes plus, two magnets. 

Electron identification (ID) was based on transverse shower shape plus match- 
ing wire chamber tracks to shower positions and energies in the electromagnetic 
calorimeter .Ej The electron ID efficiency varied from 62% below 9 GeV to 45% 
above 20 GeV. The probability to mis- ID a pion as an electron was about 0.8%. 

Muon identification was obtained from two planes of scintillation counters. The 
first plane (5.5 m x 3.0 m) of 15 counters measured the horizontal x axis while 
the second plane (3.0 m x 2.2 m) of 16 counters measured the vertical y axis. The 
counters had about 15 interaction lengths of shielding. Candidate muon tracks were 
required to pass cuts that were set using _D+ —^K fi^v decays from E791's dataJlj 
Timing from the y counters improved x position resolution. Counter efficiencies were 
measured using muons originating from the primary beam dump, and were found 
to be (99± 1)% for the y counters and (69±3)% for the x counters. The probability 
for misidentifying a pion as a muon decreased with increasing momentum; from 
about 6% at 8 GeV/c to fl.3 ± 0.1)% above 20 GeV/c. 

After reconstructionJlZI events with evidence of well-separated production (pri- 
mary) and decay (secondary) vertices were selected to separate charm candidates 
from background. Secondary vertices had to be separated from the primary vertex 
by greater than 20 a^ for £)+ decays and greater than 12 cr^ for D° and Df decays, 
where a^ was the calculated resolution of the measured longitudinal separation. In 



addition, the secondary vertex had to be separated from the closest material in the 
target foils by greater than 5 a' , where a' was the uncertainty in this separation. 
The sum of the vector momenta of the tracks from the secondary vertex was required 
to pass within 40 /ini of the primary vertex in the plane perpendicular to the beam. 
Finally, the net momentum of the charm candidate transverse to the line connect- 
ing the production and decay vertices had to be less than 300, 250, and 200 MeV/c 
for D°, D+, and D"*" candidates, respectively. The decay track candidates passed 
approximately 10 times closer to the secondary vertex than to the orimary vertex. 
These selection criteria and, where possible, the kaon identification^ requirements, 
were the same for the search mode and for its normalization signal. 

For this study a "blind" analysis technique was used. Before the selection criteria 
were finalized, all events having masses within a mass window AMg around the 
mass of D"^, D^, or D^ were "masked" so that the presence or absence of any 
potential signal candidates would not bias the choice of selection criteria. All criteria 
were then chosen by studying signal events generated by a Monte Carlo simulation 
program (see below) and background events, outside the signal windows, from real 
data. Events within the signal windows were unmasked only after this optimization. 
Background events were chosen from a mass window AM^ above and below the 
signal window AM5. The criteria were chosen to maximize the ratio Ns /^/ Nb, 
where Ns and Nb are the numbers of signal and background events, respectively. 
Asymmetric windows for the decay modes containing electrons were used to allow 
for the bremsstrahlung low-energy tail. The signal windows were: 

• 1.84 < M{D+) < 1.90 GeV/c2 for D+ -> h^fx, 

• 1.78 < M{D+) < 1.90 GeV/c^ for D+ -^ hee and h^ie, 

• 1.95 < M(D+) < 1.99 GeV/c^ for Df -^ hfj^n, 

• 1.91 < M{Df) < 1.99 GeV/c^ for Df -^ hee and Ve, 

• 1.83 < M{D") < 1.90 GeV/c^ for D" -> ^/i, 

• 1.76 < M{D") < 1.90 GeV/c2 for D^ -> ee and ^e. 

The sensitivity of the search was normalized to topologically similar Cabibbo- 
favored decays. For D"*", 24010 ± 166 D'^ — > K~tt^tt~^ decays were used; for Df, 
782±30 L»+ -^ <j>Tr+ decays were used; and for D°, 25210± 179 £>" ^ K-tt+ decays 
were used. The normalization modes widths were 10.5 MeV/c^ for Z3+, 9.5 MeV/c^ 
for D+, and 12 MeV/c^ for D^ . The upper limit for each branching fraction Bx 
was calculated using the following formula: Bx = (-/Vx/A^Norm) • (sNorm/sx) ■^Norm, 
where Nx was the 90% CL upper limit on the number of decays for rare or forbidden 
decay mode X, and ex was that mode's detection efficiency. A^Norm was the fitted 
number of normalization mode decays; eNorm was the normalization mode detection 
efficiency; and -BNorm was the normalization mode branching fraction obtained from 
the Particle Data Group.El 

The ratio of detection efficiencies is given by ENorm/^x = A^notiiiZ-^x^' where 
-^NOTin ^'^'i A^x^ ^''6 the fractions of Monte Carlo events that are reconstructed 
and pass final cuts, for the normalization and decay modes, respectively. We use 



Pythia/ jETSETiij as the physics generator and model effects of resolution, geome- 
try, magnetic fields, multiple scattering, interactions in detector material, detector 
efficiencies, and analysis cuts. The efficiencies for the normalization modes varied 
from about 0.5% to 2% and the search modes varied from about 0.1% to 2%. 

Monte Carlo studies showed that the experiment's acceptances are nearly uni- 
form across the Dalitz plots, except that the dilepton identification efficiencies typ- 
ically dropped to near zero at the dilepton mass threshold. While the loss in effi- 
ciency varied channel by channel, the efficiency typically reached its full value at 
masses only a few hundred MeV/c^ above the dilepton mass threshold. A constant 
weak-decay matrix element was used when calculating the overall detection effi- 
ciencies. Two exceptions to the use of the Monte Carlo simulations in determining 
relative efficiencies were made: those for Cerenkov identification when the number 
of kaons in the signal and normalization modes were different, and those for the 
muon identification. These efficiencies were determined from data. 

3. E791 Results 

The 90%.,CL upper limits Nx were calculated using the method of Feldman and 
CousinsE2l to account for background, and then corrected for systematic errors by 
the method of Cousins and Highland. EJ In these methods, the numbers of signal 
events are determined by simple counting, not by a fit. All results are listed in 
Table 1 and shown in Fig. 1. The kinematic criteria and removal of reflections (see 
below) were different for the D~^, D+, and D°. 

The upper limits were determined by both the number of candidate events and 
the expected number of background events within the signal region. Background 
sources that were not removed by the selection criteria discussed earlier included 
decays in which hadrons (from real, fuUy-hadronic decay vertices) were misidenti- 
fied as leptons. In the case where kaons were misidentified as leptons, candidates 
had effective masses which lie outside the signal windows. Most of these origi- 
nated from the Cabibbo-favored modes D^ — > K~Tr~^Tr~^, Df -^ K~K~^Tr~^, and 
D'^ -^ K^TT^ (and charge conjugates). These Cabibbo-favored reffections were ex- 
plicitly removed prior to the selection-criteria optimization. There remained two 
sources of background in the data: hadronic decays with pions misidentified as lep- 
tons (iVMisio) and "combinatoric" background (A^cmb) arising primarily from false 
vertices and partially reconstructed charm decays. After selection criteria were ap- 
plied and the signal windows opened, the number of events within the window was 
Nohs = NsifT + Nuism + A^Cmb- 

The background A^MisiD arose mainly from singly-Cabibbo-suppressed (SCS) 
modes. These misidentified leptons can come from hadronic showers reaching the 
muon counter, decays- in-flight, and random overlaps of tracks from otherwise sep- 
arate decays ("accidental" sources). A limit for D^ -^ K^£^£^ modes was not 
attempted, as they had relatively large feedthrough signals from copious Cabibbo- 
favored K^n^TT^ decays. Instead, the observed signals in K^l'^i^ channels were 



used to measure three dilepton misidentifieation rates under the assumption that 
the observed signals arise entirely from lepton misidentifieation. 
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Fig, 1. Final event samples for D"^ (rows 1-3), Df (rows 4-7), and D^ (row 8) decays. 
Solid curves represent estimated background; dotted curves represent signal shape of a 
number of events equal to a 90% CL upper limit. Dashed vertical lines are AMs windows. 



The following misidentifieation rates were obtained: r^^ — (7,3 ± 2,0) x 10 **, 
r^e = (2,9 ± 1,3) X 10"'', and r^e == (3,4 ± 1,4) x 10"'*, Using these rates the 
numbers of misidentified candidates, N^[^iy) (for D^ and Df) and A^misID (^O'" 
1)°), in the signal windows were estimated as follows: N^^j-^ — ^u ■ -^scs ^^'^ 
MisiD — '''u ■ ^scs' where Ng^^ and Ng^^ were the numbers of SCS hadronic 



decay candidates within the signal windows. For modes in which two possible pion 
combinations can contribute, e.g., -D"*" — » 7r+/i*/i^, twice the above rate was used. 
To estimate the combinatoric background A^cmb within a signal window AM5, 
events having masses within an adjacent background mass window AA/^ were 
counted, and this number (TVaa/b) '^^s scaled by the relative sizes of these win- 
dows: iVcmb — {^Ms/^Mb) ■ Namb- To be conservative in calculating the 90% 
confidence level upper limits, combinatoric backgrounds were taken to be zero when 
no events are located above the mass windows. In Table 1 the numbers of combina- 
toric background, misidentification background, and observed events are presented 
for all 24 modes. Also previously published results are given for comparison. 

Table 1. E791 90% confidence level (CL) branching fractions (BF) compared to previous limits. 
The background and candidate events correspond to the signal region only. 



Mode 



(Est. BG) Cand. Sys. 90% CL E791 Previous 

A'^Cmb iVMisiD Obs. Err. Num. BF Limit BF Limit 



Ref 



1.5 X 10^5 1.8 X lO"-"^ 
5.2 X 10^5 6.6 X 10^5 



D+- 
D+- 
D+- 
D+- 
D+- 
D+- 
D+- 
D+- 
D+- 



>+^ 



TT^ fi^e^ 



+ 



TT II fl 



K+ii+fi- 

K+e+e- 



1.20 
0.00 
0.00 
0.80 
0.00 
0.00 
2.20 
0.00 
0.00 



1.47 
0.90 
0.78 
0.73 
0.45 
0.39 
0.20 
0.09 
0.08 



2 
1 
1 
1 
2 
1 
3 
4 
1 



10% 
12% 
11% 
9% 
12% 
11% 
8% 
11% 
9% 



3.35 
3.53 
3.64 
2.92 
5.60 
4.05 
5.07 
8.72 
4.34 



3.4 X 10-5 

1.7 X 10-5 
9.6 X 10-5 
5.0 X 10-5 
4.4 X 10-5 
2.0 X 10-4 

6.8 X 10-5 



1.2 X 10-'' 
8.7 X 10-5 
1.1 X 10-4 
1.1 X 10-4 
9.7 X 10-5 
2.0 X 10-4 

1.3 X 10-4 



1.4 X 10-4 5.9 X 10-4 

1.6 X 10-3 
6.3 X 10-4 

1.8 X 10-4 5.9 X 10-4 

6.3 X 10-4 

6.8 X 10-4 

1.4 X 10-4 4.3 X 10-4 

2.7 X 10-4 

6.1 X 10-4 

8.2 X 10-5 4.3 X 10-4 

6.9 X 10-4 

7.3 X 10-4 
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-* K+ii+fi 


0.67 
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1.83 0.63 2 6% 3.51 5.2 x 10-^ 4.1 x IQ-^ 
1.75 0.29 9% 1.26 6.2 x 10-^ 8.2 x 10"^ 
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The sources of systematic errors in this analysis included: statistical errors from 
the fit to the normalization sample iVNorm; statistical errors on the numbers of 
Monte Carlo generated events for both N^^^ and N^'"] uncertainties in the cal- 
culation of misidentification background; and uncertainties in the relative efficiency 



for each mode, including lepton and kaon tagging efficiencies. These tagging ef- 
ficiency uncertainties included: f ) the muon counter efficiencies from both Monte 
Carlo simulation and hardware performance; 2) kaon Cerenkov identification effi- 
ciency due to differences in kinematics and modeling between data and Monte Carlo 
simulated events; and 3) the fraction of signal events (based on simulations) that 
would remain outside the signal window due to bremsstrahlung tails. The larger 
systematic errors for the Df modes, compared to the D~^ and D" modes, are due 
to the uncertainty in the branching fraction for the Df normalization mode. The 
sums, taken in quadrature, of these systematic errors are listed in Table 1. 

4. Summary 

As demonstrated elsewhere,0 one can estimate what new mass regions these results 
can probe. If one, for simplicity, assumes that gxo,X' = gw then for FCNC, using: 



mxo '^ mw 



BR{D+ -> K°fi+i^^) 
BR{D+ -> n+fi+^i-) 



(2.2) 



1/4 



(1) 



one estimates that mxo > 800 GeV/c^. Applying a similar technique for LFV, and 
now using the formula: 



mx' ~ raw 



BR{D+ ^ i?V+^^) Y^' 

BR{D+ -^ 7r+^±eT) ^ " ''J ' 



(2) 



one estimates that mx' > 650 GeV/c^, consistent with theoretical predictionsE3 

In summary, a "blind" analysis of data from Fermilab experiment E79I was 
used to obtain upper limits on the dilepton branching fractions for flavor-changing 
neutral current, lepton-number violating, and lepton-family violating decays of Z?"*", 
Df, and D'^ mesons. No evidence for any of these decays was found. Therefore, 
upper limits on the branching fractions at the 90% confidence level were presented. 
These limits represent significant improvements over previously published results. 
Eight new 13+ search modes were reported. 

For the future, FOCUS (Fermilab experiment E83I)E^ will shortly be publishing 
their results, which should improve on the E791 results. Also, because of interest 
expressedcJ in resonant decays of D^s, Fermilab experiment E79I will shortly be 
publishing results of a search for Z)° -^ V£^£^ decays. 
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